Neutron-skin thickness from the study of the anti-analog giant dipole resonance 
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The 7-decay of the anti-analog of the giant dipole resonance (AGDR) has been measured to the 
isobaric analog state excited in the p( 124 Sn,n) reaction at a beam energy of 600 MeV/nucleon. 
The energy of the transition was also calculated with state-of-the-art self-consistent random-phase 
approximation (RPA) and turned out to be very sensitive to the neutron-skin thickness (AR pn ). 
By comparing the theoretical results with the measured one, the AR pn value for 124 Sn was deduced 
to be 0.175 ± 0.048 fm, which agrees well with the previous results. The energy of the AGDR 
measured previously for 208 Pb was also used to determine the AR pn for 208 Pb. In this way a very 
precise AR pn — 0.181 ± 0.031 neutron-skin thickness has been obtained for 208 Pb. The present 
method offers new possibilities for measuring the neutron-skin thicknesses of very exotic isotopes. 

PACS numbers: 24.30.Cz, 21.10.Gv, 25.55.Kr, 27.60.+j 



Recent progress in development of radioactive beams 
has made it possible to study the structure of nuclei far 
from stability. An important issue is the size of the neu- 
tron skin of unstable neutron-rich nuclei, because this 
feature may provide fundamental nuclear structure in- 
formation. There is a renewed interest in measuring pre- 
cisely the thickness of the neutron skin, because it con- 
strains the symmetry-energy term of the nuclear equa- 
tion of state. The precise knowledge of the symmetry 
energy is essential not only for describing the structure 
of neutron-rich nuclei, but also for describing the prop- 
erties of the neutron-rich matter in nuclear astrophysics. 

The symmetry energy determines to a large extent, 
through the Equation of State (EoS), the proton frac- 
tion of neutron stars [H , the neutron skin in heavy nuclei 
and enters as input in the analysis of heavy-ion reac- 
tions [UIH, etc. Furnstahl Q demonstrated that in heavy 
nuclei there exists an almost linear empirical correlation 
between the neutron-skin thickness and theoretical pre- 
dictions for the symmetry energy of the EoS in terms 
of various mean-field approaches. This observation has 



contributed to a renewed interest in an accurate determi- 
nation of the neutron-skin thickness in neutron-rich nu- 
clei In this work, we are suggesting a new precise 
method for measuring the neutron-skin thickness using 
both stable and radioactive beams. 

In our previous work on inelastic alpha scattering, exci- 
tation of the isovector giant dipole resonance was used to 
extract the neutron-skin thickness of nuclei 

[1GI3. The 

cross section of this process depends strongly on AR pn . 
Another tool used earlier for studying the neutron-skin 
thickness, is the excitation of the isovector spin gi- 
ant dipole resonance (IVSGDR). The L=l strength of 
the IVSGDR is sensitive to the neutron-skin thickness 
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Vretenar et al. [l3j suggested another new method for 
determining the AR pn by measuring the energy of the 
GTR. Constraints on the nuclear symmetry energy and 
neutron skin were also obtained recently from studies of 
the strength of the pygmy dipole resonance [3] ■ In Ref. 
[Hj|, more nuclei were added and the theory was better 
constrained. 
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The aim of the present work is to study the energy 
of the anti-analog of the giant dipole resonance (AGDR) 
[l||, which depends strongly on the neutron-skin thick- 
ness. The non-energy-weighted sum rule (NEWSR) we 
used earlier [111, |l2[ is valid (apart from a factor of 3) also 
for the giant dipole resonance excited in charge-exchange 
reactions and predicts increasing strengths as a func- 
tion of the neutron-skin thickness. Auerbach et al. [ItJ 
derived an energy-weighted sum rule (EWSR) also for 
the dipole strengths excited in charge-exchange reactions. 
The corresponding energies are measured with respect to 
the RPA g.s. energy (IAS state) in the parent. The re- 
sult of such EWSR is almost independent of the neutron- 
skin thickness }l7j so the mean energy of the dipole 
strength should decrease with increasing dipole strengths 
and therefore with increasing neutron-skin thickness in 
consequence of NEWSR. The strong sensitivity of the 
AGDR energy on AR pn is mentioned also by Krmpotic 
18], who preformed calculations with random-phase ap- 
proximation (RPA) . In the present work, we want to use 

to 

124 Sn. 

Due to the isovector nature of the (p,n) reaction, the 
strength of the El excitation is divided into To-1, To and 
To+1 components, where To is the ground state isospin 
of the initial nucleus, which is 12 for 124 Sn. Because of 
the relevant Clebsch-Gordan coefficients [ijjj], the To-1 
component (AGDR) is favored compared to the To and 
T +l one by about a factor of T , and 2T 2 ,, respectively. 

Dipole resonances were excited earlier in (v,n) reac- 
tions at E p = 45 MeV by Sterrenburg et al. [lj| in 92 Zr, 
93 Nb, 94 Mo, 120 Sn and' in 208 Pb. Nishihara et al. 
measured also the dipole strength distributions at E p — 



such sensitivity of the energy of the AGDR on AR pn 
constrain the A.R pn of 



41 MeV. However, it was shown experimentally |2l[ |22| 
that the observed AL= 1 resonance was in general a su- 
perposition of all possible spin-flip dipole (SDR) modes 
and the non-spin-flip dipole GDR. According to the work 
of Osterfeld [l9| the spin non- flip/spin- flip ratio is favored 
at low bombarding energy (below 50 MeV) and also at 
very high bombarding energies (above 600 MeV). 

The experiments, aiming at studying the neutron-skin 
thickness of 124 Sn, were performed at GSI using 600 
MeV/nucleon 124 Sn relativistic heavy-ion beams on 2 and 
5 mm thick CH2 and 2 mm thick C targets. This allowed 
us to subtract the contribution of the C to the yield mea- 
sured from the CH2 target during the analysis. (Doing 
the experiment at 600 MeV/nucleon, in iverse kinematics 
we could increase the target thickness by a factor of 40 
compared to 50 MeV/nucleon case without loosing the 
energy resolution.) 

According to the previous experimental studies [l|l [2(| 
the excitation energy of the AGDR is expected to be at 
E x = 26 MeV. The differential cross section of the AGDR 
excited in (p,n) reaction was calculated at E p = 600 MeV 
with the computer code DW81 wave func- 

tions used by the code were constructed using the normal- 



mode formalism [25j with the code NORMOD [26j. The 
optical model parameters were taken from Ref. [271 ] . Ac- 
cording to such calculations the dipole cross section peaks 
at Q C M = 3°. 

The E x = 26 MeV and Ocm — 3° correspond in in- 
verse kinematics in the laboratory system to a scattered 
neutron with energy of E n ps 2.4 MeV and Qlab ~ 68°, 
respectively. 

The ejected neutrons were detected by a low-energy 
neutron-array (LENA) ToF spectrometer [28| , which was 
developed in Debrecen and which was placed at 1 m from 
the target and covered a laboratory scattering-angle re- 
gion of 65° < Qlab < 75°. Similar neutron spectrom- 
eters have been built also by Beyer et al. [29( and by 
Perdikakis et al. [Io| and one of those was used recently 
as an effective tool for studying Gamow^Teller giant res- 
onances in radioactive nuclear beams [3l| . 

The energy of de-exciting 7-transitions was measured 
by six large cylindrical (3.5" x 8") state-of-the-art LaBr3 
7-detectors placed at 31 cm from the target and &lab= 
21° in order to use the advantage of the large Lorentz- 
boost. The large Doppler shift (E^/Eq — 2.33) was taken 
into account in the analysis. The precise energy and ef- 
ficiency calibrations of the detectors were performed af- 
ter the experiments by using different radioactive sources 
and (p/y) reactions on different targets 32, 33|. The re- 
sponse function of the detector was also checked up to 
17.6 MeV and could be reproduced well with GEANT 
Monte-Carlo simulations. In order to make a correct en- 
ergy calibration for the AGDR, the simulations were ex- 
tended up to 40 MeV and convoluted with a Gaussian 
function with the width of the resonance. This convolu- 
tion caused about 10% lowering of the positions of the 
peaks, which was then taken into account in the cal- 
ibration of the detectors. The 7-ray energy spectrum 
measured in coincidence with the low-energy neutrons is 
shown in Fig. 1. 

The width of the peak has a value of T «3.6 MeV, 
which is much larger than the width of the previously 
measured AGDR resonances 3, 2(| ■ It can be explained 
by the Doppler broadening caused by the large solid an- 
gle of the detectors. The energy distribution of the 7- 
rays was fitted by a Lorentzian curve and a second order 
polynomial background, and the obtained parameters are 
shown in the figure. The contribution of the statistical 
error in the uncertainty of the position of the peak is 
0.2 MeV, while the systematical error coming from the 
uncertainty of the energy calibration is about 0.25 MeV 
(2.5%), which can be improved in the future. If we take 
into account the dependence of the 7-transition prob- 
ability, then the E AGDR - E IAS = 10.70 ± 0.32 MeV. 

The direct 7-branching ratio of the AGDR to the IAS 
is expected to be similar to that of the GDR to the g.s. 
in the parent nucleus, which can be calculated from the 
parameters of the GDR [l(J ■ The branching ratio ob- 
tained for 124 Sn is about 1%. In contrast, in the inves- 
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FIG. 1. The 7-ray energy spectrum measured in coincidence 
with the low-energy neutrons that fulfill the conditions of 
1.0 < E n < 3.5 MeV and 67° < Blab < 70°. The cal- 
ibrated energy scale was corrected already for the Doppler 
effect. The solid line shows the result of the fit described in 
the text. 



tigation of the electromagnetic decay properties of the 
SDR by Rodin and Dieperink 34| the 7-decay branching 
ratio was in the range of 10~ 4 . This means that in the 
n-7 coincidence spectrum the contribution of the SDR is 
suppressed by about a factor of 100. Therefore, the co- 
incidence measurements deliver a precise energy for the 
AGDR, which in case of 124 Sn a gree s well with the results 
obtained by Sterrenburg et al. [161 ] . 

The theoretical analysis employed in this work was car- 
ried out with the fully self-consistent relativistic proton- 
neutron quasiparticle random-phase approximation (pn- 
RQRPA) based on the Relativistic Hartree-Bogoliubov 
model (RHB) (H- The RQRPA was formulated in 
the canonical single-nucleon basis of the RHB model 
in Ref. [3(| and extended to the description of charge- 
exchange excitations (pn- RQRPA) in Ref. [37} . The 
RHB + pn- RQRPA model is fully self-consistent: in the 
particle-hole channel, effective Lagrangians with density- 
dependent meson- nucleon couplings are employed, and 
pairing correlations are described by the pairing part of 
the finite-range Gogny interaction [381 ] - 

For the purpose of the present study, we employ a fam- 
ily of density-dependent meson-exchange (DD-ME) in- 
teractions, for which the constraint on the symmetry en- 
ergy at saturation density has been systematically varied, 
04 = 30, 32, 34, 36 and 38 MeV, and the model parame- 



ters are adjusted to accurately reproduce nuclear matter 
properties (the binding energy, the saturation density, 
the compression modulus) and the binding energies and 
charge radii of a standard set of spherical nuclei [39} . 
These effective interactions were used to provide a micro- 
scopic estimate of the nuclear matter compressibility and 
symmetry energy in relativistic mean-field models [39} 
and in Ref. [14j to study a possible correlation between 
the observed pygmy dipole strength (PDS) in 130 > 132 Sn 
and the corresponding values for the neutron-skin thick- 
ness. In addition to the set of effective interactions with 
K nm = 250 MeV (this value reproduces the excitation 
energies of giant monopole resonances), and 04 = 30, 32, 
34, 36 and 38 MeV, the relativistic functional DD-ME2 
[40j will be used here to calculate the excitation energies 
of the AGDR with respect to the IAS, as a function of 
the neutron skin. Important for the present analysis is 
the fact that the relativistic RPA with the DD-ME2 ef- 
fective interaction predicts the dipole polarizability for 
208 Pb, an— 20.8 fm 3 , in very good agreement with the 
recently measured value: an — (20.1 ± 0.6) fm 3 

The results of the calculations for 124 Sn are shown 
in Fig. 2. The difference in the excitation energy of 
the AGDR and the IAS, calculated with the pn- RQRPA 
based on the RHB self-consistent solution for the ground- 
state of the target nucleus, is plotted as a function of 
the corresponding RHB prediction for the neutron-skin 
thickness. For the excitation energy of the AGDR we 
take the centroid of the theoretical strength distribution, 
calculated in the energy interval above the IAS that cor- 
responds to the measured spectrum of 7-ray energies: 6 
to 14.8 MeV (cf. Fig. 2). A single peak is calculated for 
the IAS. For effective interactions with increasing value 
of the symmetry energy at saturation 04 = 30, 32, 34, 
36 and 38 MeV (and correspondingly the slope of the 
symmetry energy at saturation [4l}). we find an almost 
perfect linear decrease of E(AGDR) — E(IAS) with the 
increase of the neutron skin AR np . The value calculated 
with DD-ME2 (a 4 = 32.3MeV) is denoted by the star 
symbol. 

The uncertainty of the theoretical predictions for the 
neutron-skin thicknesses is estimated to be 10 %. Such 
an uncertainty was used earlier for the differences be- 
tween the neutron and proton radii for the nuclei 116 Sn, 
124 Sn, and 208 Pb in adjusting the parameters of the ef- 
fective interactions [39|, [4(| • These effective interactions 
were also used to calculate the electric dipole polariz- 
ability and neutron-skin thickness of 208 Pb, 132 Sn and 
48 Ca, in comparison to the predictions of more than 40 
non-relativistic and relativistic mean-field effective inter- 
From the results presented in that work one 
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actions 

can also assess the accuracy of the present calculations. 

In comparison to the experimental result for 
E(AGDR) - E(IAS) we deduce the value of the neu- 



tron skin thickness in 124 Sn: AJL 



0.175 ± 0.048 fm 



(including the 10% theoretical uncertainty). In Table I, 
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this value is compared to previous results obtained with 
a variety of experimental methods. Very good agreement 
has been obtained with the previous data, which supports 
the reliability of our method. 
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FIG. 2. The difference in the excitation energy of the AGDR 
and the IAS for the target nucleus 124 Sn, calculated with the 
pn-RQRPA using five relativistic effective interactions charac- 
terized by the symmetry energy at saturation en = 30, 32, 34, 
36 and 38 MeV (squares), and the interaction DD-ME2 (ei4 = 
32.3 MeV) (star). The theoretical values E{AGDR)-E(IAS) 
are plotted as a function of the corresponding ground-state 
neutron-skin thickness AR pn , and compared to the experi- 
mental value E(AGDR) - E(IAS) = 10.70 ± 0.32 MeV. 



Fig. 3 displays the corresponding theoretical results for 
E(AGDR) - E{IAS) for the target nucleus 208 Pb, as a 
function of the corresponding ground-state neutron-skin 
thickness AR pn , and compared with the experimental re- 
sult of Sterrenburg et al. [Iff. In this case the E(AGDR) 
is calculated as the centroid of the theoretical strength 
distribution in the interval between 5 and 15 MeV above 
the IAS. The deduced value of the neutron-skin thickness 
is compared with previous results in Table II. 
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FIG. 3. Same as described in the caption to Fig. 2 but for 
the target nucleus 208 Pb. 



TABLE I. Neutron-skin thicknesses {AR pn ) of 124 Sn deter- 
mined in the present work compared to previously measured 
values. 



Method 


Ref. 


Date 


AR pn (fm) 


(P,p) 0.8 GeV 


[42,43] 


1979 


0.25 ± 0.05 


{a, a') GDR 120 MeV 


[10] 


1994 


0.21 ± 0.11 


( 3 He,£) SDR+GDR 177 MeV [12] 


2004 


0.27 ± 0.07 


antiproton absorption 


M 


2001 


0.19 ± 0.02 


pygmy res. 


[14] 


2007 


0.24 ± 0.04 


(p,p) 295 MeV 


M 


2008 


0.185 ± 0.017 


AGDR 


pres. res. 


2012 


0.175 ± 0.048 



As the 71-7 coincidence method delivered similar re- 
sults in the case of 124 Sn for the energy of the AGDR as 
obtained by Sterrenburg et al. [IH by a ToF method, it 
is reasonable to assume that their result for the energy 
of the AGDR obtained for 208 Pb is correct. 



TABLE II. Neutron-skin thicknesses of 08 Pb determined in 
the present work compared to previously measured values. 



Method 


Ref. 


Date 


AR, 




(fm) 


(p,p) 0.8 GeV 


[45] 


1980 


0.14 


± 


0.04 


(p,p) 0.65 GeV 


[46] 


1994 


0.20 


± 


0.04 


(a, a') GDR 120 MeV 


[10] 


1994 


0.19 


± 


0.09 


antiproton absorption 


m 


2001 


0.18 


± 


0.03 


(a, a') GDR 200 MeV 


[12] 


2003 


0.12 


± 


0.07 


pygmy res. 


M 


2007 


0.180 


± 


0.035 


pygmy res. 


[15] 


2010 


0.194 


± 


0.024 


(p,p") 


m 


2011 


0.156 


± 


0.025 


parity viol. (e,e) 


m 


2012 


0.33 


± 


0.17 


AGDR 


pres. res. 


2012 


0.181 


± 


0.031 



In conclusion, we have investigated the energy of the 
AGDR excited in the 124 Sn(p,n) reaction performed in 
inverse kinematics. Using the experimental results from 
this study, Ref. for 208 Pb, and RHB+pn-RQRPA 
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model, we deduce the following values of the neutron 
skin: Ai? p „=(0.175 ± 0.048) fm in 124 Sn and 0.181 ± 
0.031 fm in 208 Pb. The agreement between the Ai? p „ 
determined using measurements of the AGDR-IAS and 
previous methods is very good in both the studied cases. 
In particular, the present study supports the results from 
very recent high-resolution study of electric dipole polar- 
izability an in 208 Pb @ , respective correlation analysis of 
a d and AR pn 4.7] , as well as the Pb Radius Experiment 
(PREX) using parity- violating elastic electron scattering 
at JLAB Q. The method we have introduced provides 
not only stringent constraint to the neutron-skin thick- 
nesses in nuclei under consideration, but it also offers new 
possibilities for measuring A.R pn in rare-isotope beams. 
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